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We recently assessed the metabo- 
lism of Synechocystis sp PCC6803 
through a constraints-based reconstruc- 
tion and analysis approach and identi- 
fied its main metabolic properties. These 
include reduced metabolic robustness, in 
contrast to a high photosynthetic robust- 
ness driving the optimal autotrophic 
metabolism. Here, we address how these 
metabolic features affect biotechnological 
capabilities of this bacterium. The search 
for growth-coupled overproducer strains 
revealed that the carbon flux re-routing, 
but not the electron flux, is significantly 
more challenging under autotrophic con- 
ditions than under mixo- or heterotro- 
phic conditions. We also found that the 
blocking of the light-driven metabolism 
was required for carbon flux re-routing 
under mixotrophic conditions. Overall, 
our analysis, which represents the first 
systematic evaluation of the biotechno- 
logical capabilities of a photosynthetic 
organism, paradoxically suggests that 
the light-driven metabolism itself and its 
unique metabolic features are the main 
bottlenecks in harnessing the biotechno- 
logical potential of Synechocystis. 

Introduction 

The development of renewable energy 
sources has received significant interest 
in recent years owing to the depletion of 
fossil fuels, ever-increasing demand for 
energy, and concerns over climate change. 
A promising source of renewable energy 
is the recycling of C0 2 into usable fuels 
and fine chemicals by photosynthetic 
organisms using solar energy. There are, 



however, increasing concerns over the 
methods currently in use for producing 
biodiesel from crops and biomass. The 
problems include high production costs 
and a reduction in the amount of land avail- 
able for growing edible crops. These issues 
highlight the need for a new generation of 
biofuel technology. 1 ' 2 Cyanobacteria pos- 
sess several properties, which make them 
promising candidates for sustainable bio- 
energy generation. They are the only pro- 
karyotes capable of carrying out oxygenic 
C0 2 -fixation photosynthesis with higher 
efficiency than vascular plants. 3 " 5 The cul- 
tivation of cyanobacteria is simple, inex- 
pensive and it does not compete directly 
with agricultural crops for land or water. 
In addition, they are a source of natural 
high-value products, such as carotenoids, 
lipids, and vitamins. 6,7 They are also ame- 
nable to genetic manipulation. 8 '' These 
features have motivated recent engineer- 
ing efforts of cyanobacteria for producing 
valuable chemicals and biofuel-like com- 
pounds from the main biosynthetic build- 
ing blocks, establishing a proof of concept 
of direct biofuel production from oxygenic 
photosynthesis. 10 " 18 However, with a few 
exceptions,"' 20 the productivity has been 
very low compared with heterotrophic 
organisms. 21 

We recently reconstructed a genome- 
scale metabolic model of the cyanobacte- 
ria Synechocystis sp PCC6803 (iJN678). 22 
The model was used to study in detail 
the photosynthetic process under differ- 
ent light and inorganic carbon conditions 
as well as under genetic perturbations. 
The systems analysis identified two main 
states of the photosynthetic apparatus: 
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Figure 1. A depiction of the central metabolism of Synechocystis. Native and non-native experimentally overproduced metabolites in Synechocystis 
are represented by black and gray squares, respectively and the metabolites analyzed in this study are indicated by red lines. The carbon partition- 
ing (in %) to sugar, lipids and terpenoid biosynthesis together with the predicted carbon flux distribution (normalized to the C0 2 uptake rate) under 
autotrophic conditions is also shown. The non-native metabolites are 1-butanol (BUT), lactate (LAC), isobutyraldehyde (IBTAL). The abbreviations for 
the native metabolites are given in Nogales et a\. 22 



A C0 2 limited state (CLS) and a light- 
limited state (LLS). In addition, it was 
shown that optimal photosynthetic per- 
formance requires high photosynthetic 
robustness, including multiple lipids, 
photosynthetic pigments and alternate 
electron flow pathways (AEF), and that 
this photosynthetic robustness comes at 
a cost of reduced metabolic robustness. 
In order to explore the impact of these 
unique metabolic features and to obtain 
a better understanding of the opportu- 
nities and bottlenecks offered by cya- 
nobacteria in biotechnology, we present 
here the first analysis of the metabolic 
engineering capabilities of Synechocystis 
using iJN678. Employing an approach 
analogous to those previously used for 
in silico-driven metabolic engineering in 
heterotrophic organisms, 23 we analyzed 
how the electron and carbon flux can be 
funneled to the overproduction of both 
native (fumarate, ethanol, sucrose, and 
H ) and non-native compounds (L-lactate 
and 1-butanol). 



The above metabolites were chosen as 
representatives of key points in the metab- 
olism and/or because they have already 
been overproduced in Synechocystis (Fig. 1) . 
Growth-coupled production designs were 
attempted since they represent a stable 
phenotype, allowing for an easy selec- 
tion of the overproducing strains. 23 Flux 
balance analysis (FBA) 24 ' 25 was used to 
predict flux values in both the wild type 
and the mutants. The search for mutants 
was performed using a randomized ver- 
sion of the strategy described in Nogales 
et al. 26 Using gene-protein-reaction asso- 
ciations, which specify via Boolean rules 
the gene product(s) catalyzing a reaction, 
a mutant was created by randomly knock- 
ing out a fixed number of genes and there- 
fore, disabling flux through the affected 
reaction (s). An FBA was performed to 
determine the maximum growth rate of 
the mutant, followed by another FBA to 
determine the maximum product rate. 
Finally, a third FBA was performed to 
determine the minimum product rate, 



while enforcing the maximal growth rate; 
thus, revealing growth-coupled designs. 
By repeating this process many times, 
a list of different knockout designs was 
obtained. The number of knockouts was 
varied between 3 and 25, and the number 
of repeats was between 5 x 10 s and 5 x 
10 7 . The knockout search was performed 
under both autotrophic and mixotrophic 
conditions, employing a light-limited state 
(LLS) where the photon uptake rate was 
fixed to 30 mmol.gDW.h" 1 and a car- 
bon-limited state (CLS) where the pho- 
ton uptake rate was fixed to 100 mmol. 
gDW.h" 1 . The HCGy uptake rate was 
set to 3.7 mmol.gDW.h" 1 in all cases and 
for mixotrophic conditions, the glucose 
uptake rate was set to 1 mmol.gDW.h" 1 . 
Heterotrophic conditions were analyzed 
by setting the photon uptake rate to zero. 
The iJN678 model includes 678 genes, 
which results in a very large search space 
and to make the search more manage- 
able, a pre-processing step for reducing the 
number of target genes was performed. By 
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Table 1. Properties of the growth-coupled overproducer designs under autotrophic conditions 



Autotrophic conditions 
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(h 1 ) 


[ \ 1 II 1 1 U 1 . 
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(4C) 
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1 0.05 0.043 0.0022 


0.878 
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0.069 


0.0057 
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0.596 
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0.878 


None 








(4C) 
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Sucrose 


0.276 


None 


0.292 


None 








(12C) 


found 


found 








Lactate 


1.185 


None 


1.171 


None 








(3 C) 


found 


found 








H2 


7.154 


12 0.041 1.744 0.0723 


24.416 


8 


0.041 


19.345 


0.8001 



LLS, light limiting state; CLS, carbon limiting state; BPCY, biomass-product coupled yield. Numbers inside parenthesis represent the wild-type growth 
rate. 



identifying essential genes, genes corre- 
sponding to blocked reactions and genes, 
which participated only in conjunctive 
gene-protein-reaction rules, the number 
of target genes was reduced to 217, 222 
and 239 for autotrophic, heterotrophic, 
and mixotrophic conditions, respectively. 
The theoretical maximum production 
rate for each metabolite was estimated by 
maximizing the production of the target 
metabolite while fixing the growth rate 
to 5% of the maximal growth achieved 
under each growth condition. 

Results 

Under autotrophic conditions, we found 
that the photosynthetic states determined 
the theoretical maximum production of 
the metabolites under consideration. In 
the LLS, the yield per carbon of molecule 
of metabolites with neutral or positive oxi- 
dation states, e.g., fumarate, sucrose, and 
lactate, was significantly higher (Fig. 2; 
Table 1). Under the CLS, the yields were 
predicted to be almost identical (^0.95 
moles of carbon of metabolite/mole of 
C0 2 ), however, the theoretical maximal 
production rate of both 1-butanol and 
ethanol increased significantly suggest- 
ing that the production of the reduced 



compounds is limited by the light avail- 
ability under the LLS. The maximum 
production of H, was much higher in the 
CLS, which is consistent with the excess of 
light and the photolytic origin of the elec- 
trons used to reduce the protons. 

The search for knockout mutants 
for overproducing fumarate, ethanol, 
1-butanol, sucrose, or L-lactate under 
autotrophic conditions was very challeng- 
ing. For fumarate, a maximal produc- 
tion rate of 0.043 mmol.gDW.h" 1 in the 
LLS and of 0.069 mmol.gDW'.h-' in the 
CLS were achieved with a deletion of the 
gene slr0018, which encodes for fumarase 
(FUM) (Table 1). These predicted pro- 
duction rates are significantly lower than 
those reported in computational studies 
with heterotrophic bacteria, 23 but simi- 
lar to the in vivo yields found for other 
overproduced metabolites in cyanobacte- 
ria. 21 While cyanobacteria have success- 
fully been engineered to produce ethanol, 
1-butanol, sucrose, and L-lactate under 
autotrophic conditions, by expressing het- 
erologous enzymes, 11,13,14 our search did not 
reveal any growth-coupled mutants over- 
producing these metabolites. However, 
mutant strains able to produce small 
amounts of metabolites at their maximum 
growth rates were identified in some cases. 



For instance, the deletion of either slr2 1 32 , 
which encodes for phosphate acetyltrans- 
ferase (PTAr), or slll299 encoding for 
acetate kinase (ACKr), resulted in a theo- 
retical maximum production of ethanol 
of 0.010 mmol.gDW^.h- 1 under the LLS 
and 0.017 mmol.gDW'.h 1 under the CLS 
(Table 1). Since, pyruvate decarboxyl- 
ase and alcohol dehydrogenase II genes 
from the obligatory ethanol producing 
Zymomonas mobilis have previously been 
introduced in cyanobacteria for overpro- 
ducing ethanol, 14,27 our results suggest that 
the production rate in these recombinant 
strains could be improved by blocking 
PTAr and/or ACKr. Taken together, the 
results obtained under autotrophic condi- 
tions and the low yields found experimen- 
tally strongly suggest that a re-routing of 
the carbon flux is more difficult to achieve 
in photosynthetic organisms than in het- 
erotrophs, such as E. coli. 2i 

To give additional support to this 
hypothesis, and to investigate whether the 
autotrophic metabolism itself and/or the 
overall metabolic network of Synechocystis 
are responsible for this phenomenon, we 
searched for overproducing mutants under 
heterotrophic conditions with glucose as 
the sole carbon and energy source. Several 
growth-coupled knockouts were identified 
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Figure 2. Production envelopes for wild-type and knockout Synechocystis strains. The production envelopes for each metabolite is shown as a func- 
tion of the biomass production rate of the wild-type Synechocystis network under heterotrophic (black lines), autotrophic LLS (dark green lines), auto- 
trophic CLS (light green lines), mixotrophic LLS (red lines) and mixotrophic CLS (blue lines), as well as the growth-coupled deletion mutants identified 
(dots). The number of growth-coupled knockouts found in each condition is shown in brackets. 



for all the metabolites analyzed. In addi- 
tion, very high yields were predicted, 
ranging from 75% of the maximum pro- 
duction rate for fumarate and sucrose to 
more than 95% for ethanol, 1-butanol, 
and lactate (Table 2; Fig. 2). In fact, the 
maximum yields for ethanol, lactate, and 
fumarate were 1.9, 1.85 and 1.14 mmol/ 
mmol of glucose, respectively, in the same 
range as those predicted in silico for E. 
coli and by using a similar number of 
knockouts. 23 These findings indicate that 
it is the light-driven metabolism, rather 
than the metabolic network itself, that 
is responsible for the lack of success in 
obtaining growth coupled mutants under 
autotrophic conditions. 

Synechocystis is able to grow mixotro- 
phically with the auto- and heterotro- 
phic metabolism occurring concurrently. 
We simulated this condition in order to 
analyze the effects of the simultaneous 
presence of glucose and light on the pro- 
duction yields. The mixotrophic metabo- 
lism behaved similarly to the heterotrophic 



metabolism and provided much more flex- 
ibility for re-routing the carbon flux, com- 
pared with the autotrophic metabolism. 
High yielding growth-coupled knockouts 
were found for all the metabolites, in 
both the LLS and the CLS (Table 3; Fig. 
2). However, these yields, ranging from 
20-50% of the maximum production 
under the LLS and from 10-40% under 
the CLS, were markedly lower than those 
found under heterotrophic conditions. 
The mixotrophic mutants were found to 
share several interesting features: First, 
many equivalent overproducing mutants 
were predicted in the two photosynthetic 
states but the excess of light in CLS led 
to a significant decrease in the number of 
overproducing mutants. Under the CLS, 
the AEF pathways are essential for growth 
due to their role in redox balancing and 
they cannot be blocked simultaneously. 22 
This could indicate that the essential- 
ity of the AEF pathways under the CLS 
limits the biotechnological potential of 
Synechocystis in this state. Second, the 



blocking of key photosynthetic reactions, 
including photosystems I (PSI) and II 
(PSII), the cytochrome b ( i (CBFC) or fer- 
redoxin NADP + oxidoreductase (FNOR), 
as well as several AEFs (leading to reduced 
photosynthetic robustness) was required 
to couple the production of the target 
metabolite to growth. Third, in most of 
the cases, the overproducing mutants 
were non-viable under autotrophic con- 
ditions and glucose was used as the sole 
carbon source. Consistently and with few 
exceptions, the theoretical maximum pro- 
duction and the growth rates of the mixo- 
trophic mutants were within the range 
predicted for the heterotrophic condition 
(Fig. 2) . Fumarate was a notable exception 
and several mutants were found to be in 
the range corresponding to the LLS and 
the CLS mixotrophic states. 

In summary, overproducing mutants 
under the mixotrophic conditions were 
obtained by avoiding the light-driven 
metabolism and by reducing the photosyn- 
thetic robustness. This could indicate that 
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Table 2. Properties of the growth-coupled overproducer designs under heterotrophic conditions 



Heterotrophic conditions (0.0743) 


Metabolite 


Maximun production rate 


Number of knockouts 


Growth rate 


Production rate 


BPCY 


(mmol.gDW'h 1 ) 




(h-') 


(mmol aDW 'h" 1 ! 

\ \ i ii 1 1 \j i . y ls v v ii / 


Fumarate (4C) 


1.462 


4 


0.045 


0.719 


0.0324 


Ethanol (2 C) 


1.901 


2 


0.035 


1.218 


0.0426 


1-Butanol (4 C) 


0.951 


4 


0.036 


0.559 


0.0199 


Sucrose (12 C) 


0.443 


5 


0.062 


0.053 


0.0033 


Lactate (3 C) 


1.867 


3 


0.023 


1.444 


0.0366 


H2 


10.926 


7 


0.019 


3.195 


0.062 



LLS, light limiting state; CLS, carbon limiting state; BPCY, biomass-product coupled yield. Numbers inside parenthesis represent the wild-type growth 
rate. 



Table 3. Properties of the growth-coupled overproducer designs under mixotrophic conditions 

Mixotrophic conditions 
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0.88 


0.0533 


2.303 


6 
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0.0533 
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10 
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0.038 
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5 


0.03 
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0.0088 


0.767 


5 


0.03 


0.295 


0.0088 
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6 
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0.0291 
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0.0291 


H2 
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11 


0.112 


2.897 


0.3248 
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8 


0.156 
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LLS, light limiting state; CLS, carbon limiting state; BPCY, biomass-product coupled yield. Numbers inside parenthesis represent the wild-type growth 
rate. 



the photosynthetic processes do indeed 
limit the possibilities in re-routing of the 
carbon flux and consequently the overpro- 
duction of the target metabolites. 

Photohydrogen production has pre- 
viously been reported in Synechocystis 
and other cyanobacteria. 1011 In order to 
explore how the electron flux can be re- 
routed under different growth conditions, 
we extended our analysis to search for H 2 
overproducing mutants. Several mutants 
with high production rates were found in 
all the growth conditions. The yields of 
H 2 were highest under auto- and mixo- 
trophic conditions and the excess of light 
under the CLS increased the production 
rates almost 4-fold (Tables 1-3). This was 
expected since the photolysis of H 2 0 is 
the main source of electrons in presence 
of glucose and it is the sole electron source 
under autotrophic conditions. An interest- 
ing finding was that the H 2 overproduc- 
ing mutants were achieved by blocking 



the AEF pathways almost exclusively, 
thus reducing the photosynthetic robust- 
ness. Consequently the H 2 production 
remained as the main electron sink in the 
network (Tables 1 and 3). 

This study is a first step toward evalu- 
ating how the particular metabolic fea- 
tures of Synechocystis that we revealed in 
our previous work, may affect the bio- 
technological potential of this organism. 
It must be noted that since the random 
search strategy is unable to test all pos- 
sible knockouts (unless the number of 
knockouts is very small), the existence 
of growth-coupled autotrophic mutants 
for some of the compounds under inves- 
tigation cannot be ruled out. The mul- 
tiple overproducing mutants found under 
heterotrophic and mixotrophic condi- 
tions indicate, however, that the lack of 
growth-coupled autotrophic mutants is 
rather associated with the properties of 
light-driven metabolism than the search 



strategy itself. Another caveat is the lim- 
ited number of compounds examined in 
this study and the conclusions may not 
apply to other metabolites or compounds 
of interest. In addition, it must be taken 
into account that we have mainly explored 
the growth-coupled biotechnological 
capabilities of Synechocystis and that alter- 
native production strategies were beyond 
the scope of this study. Nevertheless, the 
analysis of the Synechocystis network under 
different growth condition presented here 
and the comparison of the results with 
those obtained previously with hetero- 
trophic bacteria, offers a general view of 
the biotechnological capabilities of this 
cyanobacterium. 

Several conclusions can be inferred 
from this study. First, we have found 
that the re-routing of the carbon flux 
in Synechocystis under autotrophic con- 
ditions is significantly more challeng- 
ing than under hetero- or mixo-trophic 
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conditions. The constrained carbon flux 
distribution, in which up to 80-85% of 
the total fixed CO, is funneled to sugar 
biosynthesis 18 together with the metabolic 
peculiarities of Synechocystis under auto- 
trophic conditions, 22 could be the pri- 
mary contributing factors. Mixotrophic 
conditions could be used in order to 
bypass this limitation, while allowing net 
C0 2 fixation. However, we found that the 
blocking of the light-driven metabolism 
and that reduction of the photosynthetic 
robustness was a prerequisite for coupling 
the production of the target metabolites 
to growth. In addition, the mixotro- 
phic metabolism in these mutant strains 
resembled the heterotrophic metabolism, 
with (net) CG* 2 fixation absent in the 
most of the cases. Finally, we have noted 
that, in contrast to the carbon flux, the 
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